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Abstract 



H ' The effect of anomalous chromomagnetic (/i) and chromoelectric couplings id) of the 
gluon to the top quark are considered in e^e^ tt, with unpolarized and longitudinally 
polarized electron beams. The total cross section, as well as t and t polarizations are 
calculated to order as in the presence of the anomalous couplings. One of the two linear 
combinations of t and t polarizations is CP even, while the other is CP odd. The limits 
that could be obtained at a typical future linear collider with an integrated luminosity of 
50 fb^^ and a total cm. energy of 500 GeV on the most sensitive CP-even combination 
of anomalous couplings are estimated as —3 < Re(/x) < 2, for Im(yu) = = d and 



Im(/i)^ + < 2.25 for Re(/i) = 0. There is an improvement by roughly a factor of 2 
at 1000 GeV. On the other hand, from the CP-odd combination, we derive the possible 



complementary bounds as —3.6 < lm{fi*d) < 3.6 for lm{d) = and —10 < lm{d) < 10 
for lm{fj,*d) — 0, for a cm. energy of 500 GeV. The corresponding limit for 1000 GeV 
is almost an order of magnitude better for Im(//*d), though somewhat worse for 
lm{d). Results for the cm. energies 500 GeV and 1000 GeV, if combined, would yield 
independent limits on the two CP-violating parameters of —0.8 < Im{ij,*d) < 0.8 and 
-11 < lm{d) < 11. 

PACS numbers: 12.38.Bx, 13.40.Em, 13.88.+e, 14.65.HA 
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1 Introduction 



The discovery of a heavy top quark, with a mass of = 175 ± 6 GeV which is far 
larger than that of all other quarks, opens up the possibility that the top quark may have 
properties very different from those of the other quarks. Observation of these properties 
might even signal new physics beyond the standard model. Several efforts in the past few 
years have gone into the investigation of the potential of different experiments to study 
possible new interactions of the top quark. In particular, possible anomalous couplings of 
the top quark to electroweak gauge bosonsQ and to gluons 0, |^ have also been discussed. 
Top polarization is especially useful in such studies |^, |§, 0, §] , because with a mass around 
175 GeV, the top quark decays before it can hadronize @, and all spin information is 
preserved in the decay distributions. 

In this paper, we investigate the potential of e^e~ experiments at a future linear collider 
with centre-of-mass (cm.) energies of 500 GeV or higher, to study anomalous chro- 
momagnetic and chromoelectric dipole couplings of the top quark to gluons. So far, a 
considerable amount of earlier work on the topic of anomalous gluon couplings has concen- 
trated on hadron colliders. But also high energy e^e~ experiments with sufficiently high 
luminosities would provide a relatively clean environment to probe the standard model 
for anomalous gluon couplings. While earlier efforts in the context of e^e~ colliders are 
mainly based on an analysis of the gluon distribution 0] in e^e~ — tig, we look at the 
possible information that could be obtained from studying the total cross section, and the 
polarization of t and t separately. This has the advantage over t and t spin correlations 
that, because the polarization of only one of t and t is analyzed by means of a definite 
decay channel, the other is free to decay into any channel. This leads to much better 
statistics compared to the case when t — t spin correlations are considered, where definite 
t and t channels have to be used as analyzers. 

We find that there are three independent quantities, viz., the cross section, and the CP- 
even and CP-odd linear combinations of the t and t polarizations, which can be used to 
^References to the voluminous literature on this subject can be found, for example, in [Q. 
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probe separately the CP-even chromomagnetic dipole moment coupling and the CP-odd 
chromoelectric dipole coupling. Of these, the CP-odd combination provides the most 
sensitive probe of the imaginary parts of the anomalous couplings. 

We obtain here the 90% confidence level (C.L.) limits that would be possible at a typical 
e'^e~ linear collider, with integrated luminosity 50 fb~^. We have considered two possible 
centre-of-mass (cm.) energies, viz., 500 GeV and 1000 GeV. We have also considered the 
effect of beam polarization on the sensitivity of the measurements. 

Our results on the constraints on the CP-violating chromoelectric dipole moments were 
presented in Ref. . We include in this paper also the limits that would be obtainable 



on CP-conserving chromomagnetic dipole moments, and combinations of couplings. 

The paper is organized as follows: In Section 2, we introduce the effective action for 
the general massive anomalous qqg vertex. The subsequent discussion presents the full 
framework necessary for obtaining the analytical extensions that modify the standard 
QCD one-loop predictions for quark-antiquark production. We then give explicit combi- 
nations of the total polarized cross sections which are sensitive to the chromomagnetic 
and chromoelectric dipole moments, respectively. Next, Section 3 focuses on the numer- 
ical estimates of these observables for various kinematical regions readily accessible at a 
future e~^e~ collider. The final Section 4 is the Conclusion. 

2 Calculation of cross section and top polarization 
with anomalous couplings 

An effective ttg vertex can be written in the form 
where 



nit fnt 



= iA"; (2) 
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A'^ being the SU{3) Gell-Mann matrices, and and are respectively the momentum 
and polarization four-vectors of the gluon. This is the most general Lorentz- and colour- 
invariant trilinear coupling (additional quadrilinear terms are needed for local colour in- 
variance, but we do not need them here). The fi and d terms are the chromomagnetic 
and chromoelectric dipole terms, respectively. These dipole couplings fi and d are in fact 
momentum-dependent form factors, and complex in general. They parameterize the ef- 
fects arising at the loop level, which in principle could arise from both, interactions within 
the Standard Model (SM) and possible new interactions. If the measured values of the 
corresponding form factors deviate from the theoretical predictions of the conventional 
SM, this would indicate the presence of additional, new physics. Note that a possible 
experimental determination of non-vanishing Im(/i) [Im((i)] does not necessarily imply 
the existence of CP-conserving [CP-violating] new physics. 

We use Eq. (P to calculate the ti total cross section and the t and t polarizations to 
order = g'^/lAn). The extra diagrams contributing to this order are shown in Fig. 1, 
where the large dots represent anomalous couplings. The anomalous couplings enter 
in the amplitudes for soft and collinear gluon emissions. The infrared divergences in 
the amplitudes for soft gluon emission and in the virtual gluon corrections cancel as in 
standard QCD, since the anomalous terms vanish in the infrared limit. Moreover, we do 
not include anomalous couplings in the loop diagrams. The latter are included merely to 
regulate the infrared divergences. We thus study how anomalous couplings at tree-level 
would modify standard QCD predictions at order as- 

For heavy-quark production, the standard QCD one-loop corrections to the total cross 
section and the longitudinal spin polarization [|lT| were calculated in closed analytic form 



before. Those results have been extended here to the case when anomalous couplings are 
present. 

The total unpolarized production rate is given only in terms of the VV and AA parity- 
parity combinations for the Born contributions: 

cTBorn (e^ , Z qq) = - t;^)^^^ + v^a^^, (3) 



where the mass parameters are v = \/l — ^ and ^ = Arrig/s. The 0{as) unpolarized case 
has the cross section 



a (e+e- ^^,Z^qq) = ^3 - v^)a^^ c""^ + v^'a^^c^^ 



(4) 



where the VV and AA factors that multiply with the appropriate Born terms are given 
below. 



.yv 



a. 



1 + —C 

2tt ' 



r - V-— In ]--h- -h + ^77— 7T + ^5 + Aq 



(5) 



In this equation, the contribution from the virtual gluon loop is denoted as 



+3t;ln (^^ ) -4. 



Li2 



2v 



l-v 



Li2 



2t; 
l + v 



(6) 



Here, the qqg phase-space integrals are abbreviated by Jj, and specify the results after 
the (soft) IR divergences have canceled. The explicit analytical expressions for these 
phase-space integrals may be found in |T^. The additional component stemming from 
the anomalous gluon bremsstrahlung is given by 

2 



^0 



(2 + 0^^ 

For the AA contribution we find: 



Re(/i)(/i + ^4) + ^ [W + W) ih - 2/8) 



(7) 



r + 2- In I I + — Ji 12 13 H 5-/4 h + An 



V VI — f/ f'^ V V 



with the following anomalous part 



AAA 

^0 



-| Re(/i){ - (4 - 0/1 + (2 + 0/4} + (l/^P + MP) I (I + ^ - /i + e/4 



4 4 ' 

--(2-0^8 + ^/9 



(9) 



The remaining VA and parts are identical and only contribute to the spin-dependent 
cross section. In the absence of anomalous couplings {fi = d = 0), the cross section for 
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longitudinally polarized quarks of helicity ±| is given by 

a (e+e- ^ 7, Z g(A±) q) = it;(3 - t;2)a^^c^^ + t;V^'^c^^ ± (10) 



The multiplication factors are expressed in terms of phase-space integrals of type Si. 

(11) 



where 



VA,AV 



(4 - 0^1 - (4 - 50^2 - 2(4 - 30^4 - e(l - 0('53 + + - Sr 



-2S8 + (2 - 0^9 + (6 - O'^io - 2^11 + 2(1 - 0(2 - OS 



12 



121 



The full analytic forms of all relevant 5* integrals are too lengthy to be exhibited here. 
Most of them are compiled in Ref. [jl2l, except for the four additional integrals 514, 5*15, 
5*16, and 5*17, which are discussed in more detail in the appendix. 

Including spins for the quark or the antiquark introduces additional spin-flip terms in the 
0{as) c factors given in Eqs. (|^), (|^) and (|ll|). For longitudinal quark polarization we 
find 



— 1 



27r 



(13) 



The individual parity-parity combinations are 



(2 + O^A^ 



vv 



8 Im (/i*c/) 



2 l-j s. 



4\ 

1 — ^ I 'S's — 5*9 — 5*10 



S 
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+2 h 5i3 + |(5i5 + 5i7) 



+ Im(c?) 



8(1 - 0^1 - {8 - 3^2 - 0}52 + {8 + e(2 - 30}5'4 
-e(2 + 0^6 - 4^8 + 4(1 - 0^9 - 4(3 + 0^10 + 4^11 
+e(2-0^i4 



(14) 



4Im (/i*(i) 



-(1-0(2-0^1 



5-- .3 



;i- 0(^9 + ^10) + ^11 



+2 U - -\ Sr, - 2 fl - I) ^15 - ^5i6 - 2 (1 - ^17 



+ Im{d) 



2^S, - (1 - 0(4 - 0iS2 - S,) - e(l - OSe - (2 - OSi 



+(2 - 30^9 - (6 - 50^10 + (2 + OiSu + 2^13) + e(l - OSu 



, (15) 



= Re(/i) 



+ ^7) - 2Ss + (2 - 0(^9 + ^10) - 2Su 



±2ilm(/i) 



(2 - - CS_ 
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2S 
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(16) 



- (4 - ^)Si - ^87 + ASs + + (4 - O-^io - 4^11 



Using charge conjugation in the final state, one can readily obtain the corresponding 
expressions for (longitudinal) antiquark polarization. In the following, we denote the 
antiquark results by an additional bar, i.e. A^: 



where the following identities hold 



^5 



^5 > 



.AA 



'5 ' 



(17) 
(18) 



with 



XAV 



.ava 
^5 • 



(19) 

(20) 



Considering the above expressions, we can construct the following combinations of polar- 
ization asymmetries of t and t, 



A<j(+) 



A<j(-) 



(21) 
(22) 



a(t)-a(i)-a(T) + a(i) 

>(T)-<^(i)+^(T)-^(i); 

where (7(f), a(t) refer respectively to the cross sections for top and antitop with positive 
hehcity, and cr{l), a{l) are the same quantities with negative helicity. Of these, Aa^'^^ is 
CP even and Aa^~^ is CP odd. This is obvious from the fact that under C, a and a get 
interchanged, while under P, the helicities of both t and t get fiipped. Consequently, a 
and Aa^^\ both nonzero in standard QCD, receive contributions from combinations of 



anomalous couplings which are CP even, viz., Im(yu)^ + and Re(/i). On the other hand, 
Acr^^) vanishes in standard QCD, and in the presence of anomalous couplings it depends 
only on the CP-odd variables lm(n*d) and lm(d). That Aa^~^ depends on the imaginary 
parts rather than the real parts of a combination of couplings follows from the fact that 
it is even under naive time reversal T^, i.e., reversal of all momenta, without change in 
helicities, and without interchange of initial and final states (as would have been required 
by genuine time reversal). As a consequence, it is odd under CPTn, and imaginary parts 
of couplings have to appear in order to avoid conflict with the CPT theorem. 



3 Numerical results 



Figs. 2a and 2b show the dependence of Acr'^+) on ylm(yu)2 + \d\'^ and Re(yu) at ^/s = 500 
GeV and 1000 GeV, respectively. Figs. 3a and 3b depict Aa^^^ plotted against lm{n*d) 
and lm{d). 

We would now like to see how sensitive experiments at a future linear collider would be to 
the anomalous quantities /i and d. Since these are determined through Ao"^^-* and Aa^~\ 
we should know how well these latter quantities can be measured experimentally. The 
quantity cr{]) — cr(J,) is simply the polarization Pf of the top quark in the production 
process, and it can be determined by looking at the angular distribution of the decay 
products of the top. For example, in the decay mode t W~^b, the angular distribution 
of the b quark in the t quark rest frame is 

1 dV 



^[l + P,ap,-s), (23) 



r dQb 47r 

where pi, is the 6-quark momentum and s is the top spin. This angular distribution can 
be used to determine Pf. In Eq. (|23D, the parameter a {\a\ < 1) is a constant known as 
the analyzing power for the decay channel. In this particular case, a is given by 

a = i ^ f . 24 
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There is an analogous expression for the decay t — W~b. If the angular distribution of a 
lepton or a jet arising from W decay is used to determine Pt, the corresponding analyzing 
power would be different. The efficiency with which the top or antitop polarization can 
be measured will depend not only on the analyzing power of the channel, but also on the 
detection efficiency of the observed particles, like the b quark. While the tagging efficiency 
of the b would be much better at e^e~ colliders than at hadron colliders, because of the 
lower hadronic activity, it would nevertheless depend on the degree to which backgrounds 
are understood and eliminated. A more detailed analysis in the context of specific ex- 
perimental conditions would be needed to obtain the overall top polarization detection 
efficiency, which we parameterize as e in what follows. 

We use our expressions to obtain simultaneous 90% confidence level (CL) limits that 
could be obtained at a future linear collider with an integrated luminosity of 50 fb^^. We 
do this by equating the magnitude of the difference between the values for a quantity 
with and without anomalous couplings to 2.15 times the statistical error expected. Thus, 
the limiting values of Im(/i)^ + \d\'^ and Re(/i) for an integrated luminosity L and a top 
detection efficiency of e are obtained from 



eL\a{i^,d) - Osm \ =2.15 \/LasM (25) 

and 



eL AaW(/x,rf)-A(TW = 2.15 ^/i^ | a,,,(T) - a,„(T) | , (26) 
whereas the limiting values of lm{fi*d) and lm{d) are obtained from 

eL I A(t(-) ifi, d) - AairJ \=2.15^L\ (t,„(T) + asM(T) I • (27) 

In the above expressions, the subscript "SM" denotes the value expected in the standard 
model, with fi = d = 0. We use L = 50 fb~^ and e = 0.1 in our numerical estimates. 
The value of e depends on the details of the detector, as well as the kinematical cuts 
employed. The value we use is only representative, and it would be easy to obtain limits 
for any other value of e by appropriate scaling, because of the simple dependence on e in 
Eqs. (H), (13), and (|2|). 
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For the running of the strong couphng, we choose a^^\Mz) = 0.118 (with Mz = 
91.178 GeV) in the modified minimal subtraction scheme and use the appropriate condi- 
tions to match for six active fiavoursQ. 



Eqs. (25) and (p6|) are used to obtain contours in the plane of y Im(/i)^ + \d\'^ and Re(/i). 
The contours from (p5D are shown in Figs. 4a and 4b for ^/s = 500 GeV and 1000 GeV, 
respectively. The contours obtained from (pB]) are shown in Figs. 5a and 5b. Eq. (p7| ) 
gives contours in the plane of lm{^*d) and lm{d), which are displayed in Figs. 6a and 
6b. For different e~ longitudinal beam polarizations P_, the corresponding contours are 
presented in Figs. 4-6. In Figs. 4 and 5, the allowed regions are the ones shown below the 
respective contours. In Figs. 6a,b, the allowed regions are bands lying between the upper 
and lower straight lines. 

In principle, the measurement of two independent quantities a and Aa^~^^ could have 
given independent limits on Im(/i)^ + and Re(/i). However, a comparison of Figs. 4 
and 5 shows that this possibility is not realized. A superposition of Figs. 4a,b and 5a,b is 
shown in Fig. 7a. It can be seen that the allowed region coming from a lies entirely within 
the allowed region from Aa^'^\ with no intersection between the two. Thus a is more 
sensitive to the anomalous CP-even couplings as compared to Aa'-^-'. Moreover, each 
contour for ^/s = 1000 GeV lies within the corresponding contour for y/s = 500 GeV, 
showing a uniform increase in sensitivity with cm. energy. 

The other conclusion that can be drawn from Figs. 4 and 5 is that a large left-handed 
polarization leads to increase in sensitivity. 

The best limits obtainable from a at -^i = 500 GeV are —3 < Re(/i) < 2 for Im(/i) = = d, 



and y Im(/i)2 + < 2.25 for Re(/i) = 0. These limits are improved by about a factor 2 
in going to -^/i = 1000 GeV for the same integrated luminosity. 

In the case of CP-odd combinations of the couplings lm{fi*d) and lm{d), there is only one 
^It is common to indicate the number of active flavours as superscript of as- For practical purposes 
one usually selects bottom production as reference. Here, our choice for af \Mz) translates to af\Mt = 
172.1GeV) = 0.10811 
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measurable CP-odd quantity Aa^~^ at each cm. energy, and therefore independent hmits 
on the two combinations are not possible. Fig. 6a shows that for y/s = 500 GeV, the best 
limits which can be obtained are for P = —1, viz., —3.6 < lm.{fi*d) < 3.6 for lm{d) = 
and — 10 < Im((i) < 10 for lm{fi*d) = 0. The corresponding limits for y/s = 1000 GeV, 
as seen from Fig. 6b, are —0.4 < Im(/i*d) < 0.4 for lm{d) = and —20 < lm{d) < 20 
for Im(/i*(i) = 0. However, if a measurement of Aa^^^ is made at two cm. energies, a 
relatively narrow allowed range can be obtained, allowing independent limits to be placed 
on both Im(/i*(i) and lm{d). This is demonstrated in Fig. 7b. In fact, the improvement 
in the limit on lm{fi*d) in going from y/s = 500 GeV to y/s = 1000 GeV is considerable. 
The possible limits are 

- 0.8 < Im(/x*d) < 0.8, -11 < lm{d) < 11. (28) 

These limits may be compared with the limits obtainable from gluon jet energy distribu- 
tion in e^e^ —>■ tig @]. While our proposal for the CP-even case seems to fare worse, for 
the CP-odd case, our proposal can be competitive. It should however be emphasized that 
in the case of the CP-odd couplings, we are proposing the measurement of a genuinely 
CP- violating quantity, whereas the analysis in 0] is merely based on the energy spectrum 
resulting from both CP-odd and CP-even couplings. In case of Aa^~\ the dependence on 
e~ beam polarization is rather mild. 

4 Conclusions 

It is worthwhile noting that we have used a rather conservative value of e = 0.1 for top 
detection and polarization analysis. A better efficiency would lead to an improvement in 
the limits, as would a higher luminosity. 

We have not considered the effect of initial-state radiation in this work. We have also 
ignored possible effects of collinear gluon emission from one of the decay products of t or 
t. A complete analysis should indeed incorporate these effects, as well as a study of t and 
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t decay distributions which can be used to measure the polarizations. However, we do not 
expect our conclusions to change drastically when these effects are taken into account. 

In summary, we have examined the capability of total cross section and single quark 
polarization in e+e~ ti to measure or put limits on anomalous chromomagnetic and 
chromoelectric dipole couplings. While the total cross section measurements can give, 
for the luminosities assumed, limits of order 1 on the CP-even couplings, only the CP- 
violating combination of top and antitop polarizations is sensitive to anomalous couplings, 
and can yield a limit of the order of 1 on a CP-odd combination of anomalous couplings. 
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Appendix: Additional Phase-Space Integrals 



The complicated three-body phase-space for the process e^e^ (l{])(ig is best solved 
analytically by using the kinematic variables y = 1 — pi ■ g/g^ and z = 1 — p2 ■ q/q"^, 
which are natural dimensionless parameters referring to the quark and antiquark energies 
(including the radiated gluon) in the centre-of-mass system. Then, the relevant phase- 
space integrands all result in simple rational functions containing polynomials in y and z, 
and the corresponding integration boundary is described by the symmetric solution of a 



{y, 2;)-biquadratic form |n|, |T2|. 



Apart from the detailed integral list of Ref. [|T2|, in this particular calculation four new 
integrals emerged. For completeness, we give here the following full analytical results: 



S 



14 



dy dz 



2 2 + y/^ 

e 2(2 - v^) 



ln(2-v/e) + |lne 



1 

2' 



(A.l) 



S 



15 



dy dz 



|lne-ln(2- V^) 



(A.2) 



S 



16 



dy dz 



y'^z 



^|^^ + ^(^e-l)ln(2-V?)-ie(l + ie) 



+ 512^ 4-e + (16-301ne 



^4 13 



If + ±(A + 1. 

2^ ^ 16^ / ^ 24 



24' 



(A.3) 



5* 



17 



dy dz 



yz 



14 



^e'(6 - [ - I + ln(2 - V^)J + ^(4 - 



(2 - 



(A.4) 



where the usual mass parameters are v = y/1 — ^ and ^ = 4m^/s. 



Note that there are no soft divergences in these particular integrals. The coUinear diver- 
gences contained in these integrals arise from the massless character of the fermion field, 
and are easily identified by observing the limit ,^ ^ 0. In our full analytical expressions 
for the polarized cross sections with VV and AA parity-parity combinations, integral 5*14 
is multiplied by the mass factor ^, viz. Eqs. ([I^) and (0). This produces additional finite 
contributions, which originate from a coUinear helicity-fiip with anomalous chromoelectric 
couplings, and are absent in a naive massless model. 
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Figure Captions 



Fig. 1: Additional Feynman diagrams contributing to a {e^e" ^, Z ^ tt) that ac- 
count for anomalous gluon couplings at 0{as). The large dots represent 
anomalous ttg insertions according to the effective action Eq. (1). 

Fig. 2: Surface plots displaying the dependence of the polarization asymmetry Aa^'^^ 



on Y'lm(/i)-^ + MP ^iid Re(yLi) with initial electron beam polarization P_ = — 1 
and cm. energies (a) = 500 GeV, (b) ^/s = 1000 GeV. 

Fig. 3: A(7^~^ surface plots showing the hnear dependence on Im(//*d) and lm{d) for 
(a) ^/s = 500 GeV and (b) ^/s — 1000 GeV with longitudinal beam polariza- 
tion P_ = -1. 

Fig. 4: Contour plots showing the allowed regions for with 90% confidence level 
(integrated luminosity L = 50 fb~^ and top detection efficiency e = 0.1). 
Representative cm. energies are (a) ^/s = 500 GeV and (b) ^/s = 1000 GeV 
for various longitudinal electron polarizations. 

Fig. 5: A(7(+^ contour plots with 90% confidence level at cm. energies: (a) a/s — 
500 GeV and (b) = 1000 GeV. 

Fig. 6: Aa^~^ contour plots with 90% confidence level at cm. energies: (a) ^/s — 
500 GeV and (b) y/s = 1000 GeV. 

Fig. 7: (a) Superposition of Figs. 4 and 5, displaying combined allowed regions for 
both and Aa*^^-' polarization asymmetries, (b) Intersecting area resulting 
from two independent Aa^~^ measurements at ^/s — 500 GeV, 1000 GeV. 





Figure 1 




Figure 2a 




Figure 2b 




Figure 3a 




Figure 3b 




Figure 4a 




Figure 4b 




Figure 5a 




Figure 5b 




Figure 6a 




Figure 6b 




-20 -15 -10 -5 5 10 15 20 

lm(d) 



Figure 7b 



